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DYNAMICS OF A RIGID BODY 



3. 1 INTRODUCTION AND DEFINITION. 



3. 1. 1 In general, the methods of particle dynamics are insufficient to adequately 
1 describe the motion of a physical body. In this section, an idealized but useful treat- 

ment of rigid-body motion is developed and extended from the concepts of particle 
T motion presented previously. The definition of a rigid body is of paramount importance; 

one must constantly bear in mind that the theorems and laws developed in this section 
1 are fundamentally dependent on this definition. 

A rigid bod y is defined as a grouping of particles of masses m. , m„ , m , . . . , m 
I - 1 A 3 n 

| such that the distance between any two of them does not change with time. Note that the 

possibility of internal forces acting within the body is not precluded. Rather, one can 
J infer from the definition that the internal forces are such that the distance between any 

two points remains invariant with respect to time. Thus, deformations, such as those 
due to stresses and strains, are not considered in this presentation. Furthermore, it 
is assumed that the rigid body can move freely in space according to the forces im- 
pressed on it; physical constraints on its motion are not considered. 

3. 2 INDEPENDENT PARAMETERS. 

3. 2. 1 We may now raise the question: How many independent parameters are required 
to specify the instantaneous position and orientation of a rigid body with respect to a set 
of coordinate axes? Obviously, knowing the location fx_ . v. . * ■* nf nno ™,.+,-„i, 
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Nor would the locations of two points P 1 and P 2 (Fig. 31) indicate the orientation of 
the body. 
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suppose the instantaneous locations of three points P p and 

J. Cj 



P„ of the 



However, 

body were known: The distances between these points (Fig. 32) are fixed by definition. 

Let d be the distance between P and P 
d 2 the distance between P and P , and 
d 3 the distance between P and P . Then 
the coordinates of P x , P g and P g , namely 




(x 3 



z l> ' < X 2 



y 2 > z 2 ) - and 



Y 3 » z 3 ) , respectively, would suffice 
to specify the position and orientation of the 
body at that time. Note, though, that the 
orientation is known only if the three points 
are not colinear. 



Fig. 32 



These nine coordinates are not independent, as we shall see; that is, the body does not 
have nine independent degrees of freedom of motion. For example, three independent 
parameters are required to specify the location of one point, say P . However, since 
the distance A % between P 1 and P g remains fixed (by definition of rigid body) , P 
is thus located on a sphere of radius d 1 about P as center. Hence, only two inde- 
pendent parameters are required to locate P , once P is specified. Similarly, the 
third point P g is on a sphere of radius d g about P and on another sphere of radius 
d 2 about P 2 . From which, we see that P is on the circle of intersection of both 
spheres, and, therefore, only one independent parameter is required to locate it after 
the locations of P 1 and P g have been determined. A maximum of six independent 
parameters are thus required to specify the location of the three points and, hence, the 
position and orientation of the rigid body to which they are attached. Six independent 
parameters are required wheu the body can move freely in space subject to forces acting 
on it. Fewer independent parameters may suffice when physical constraints are placed 
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One may also conceive of these six parameters as being divided Into two groups: three 
parameters are required to specify the location of one point in the body with respect to 
some set of axes, and three parameters are required to specify the orientation of the 
body with respect to that point. The body is then said to have six degrees of freedom 
of motion. 
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In Fig. 33, X^ , Y^ , and Z^ are the axes of an inertial coordinate system with origin 
Q . Let X Q , Y Q r and Z Q be another Cartesian coordinate system whose axes always 
remain parallel to the respective axes X^ , Y^ , and Z' , but whose origin O is a 
point of the rigid body. The X Q , Y Q and Z Q axes will be called the translating refe r- 
®S2£JZ£tem. Let X , Y , and Z be a Cartesian coordinate system rigidly attached 
to the body, whose origin O is the same as that of the translating reference system; 
these axes will be known as the body axes , and, of course, will rotate with the body. 
The configuration of the rigid body B can then be described with respect to the body 
axes X , Y , and Z. 
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The motion of the body is composed of: 

• The translation of the origin O of the body axes 

• The rotation of the body axes X , Y , and Z about the translating 
reference frame /X , Y , Z \ 

Thus, it is necessary to find a method which gives the instantaneous orientation of the 
X , Y , and Z axes with respect to the X Q , Y q , and Z q axes. For the moment, 
however, we turn to the translation problem; that of specifying the motion of the origin 
of the body axes. 

3. 3 LINEAR MOTION. 

3.3.1 Returning to the original concept of a rigid body (i. e. , an assemblage of particles 

of mass m, , m , m c ,... , m ) , let r. , r , and f be the corresponding radius 
i o o n l & o 

vectors of the points from the origin of some inertial coordinate system. Let ( F ext ) j 
and ( F. .). represent the sums of all external and internal forces acting on the 
i particle, where the internal forces are attributed to attractions and repulsions of 
other particles. Then, Newton's second law gives 

m A = < ¥ ext>i + < ¥ int>i i3 - 1} 

Summing over all the particles of the body, 



f m i ? i = f ( F ext>i + !< F int>i < 3 - 2) 



Since the internal forces are due to attractions and repulsions between the particles and 
are not the result of external causes, Newton's third law holds; that is, if two particles 
exert forces on each other, the force exerted by the first on the second is equal to and 
opposite to the force exerted by the second on the first. The sum of all the internal 
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It is now convenient to define a point in the body known as the center of mass. The 
position vector r of the center of mass with respect to some arbitrary coordinate 
system is given by 



cm 



cm 



cm 



2 m.f. 
i 1 l 


2 m.r. 
i l 1 
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2 m.x. 
i i i 

2 m.y. 
i r i 

2 m.z. 

_i 1 * 


2 m. 
i i 



(3.4) 



where m is the total mass 2 m . Note that if the origin of the coordinate system is 
taken to be the center of mass itself, then the coordinates are (0,0,0). Therefore, 



2 m.r. = 
i i i 



(3.5) 
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Differentiating Eq. (3.4) twice with respect to time, we obtain 



* =5^f»iV 



and, since each particle is considered to have constant mass m. , 

l 



2 m.r. 



l l 



r = 



m 



Substituting from Eq. (3.3) and multiplying by m , we have 
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mr = 2 (F ,). 
* x ext'i 



(3.6) 
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Thus, the mass center of a system of particles moves as if all the mass were concen- 
trated at that point and as if all external forces acted there. Hence, the component 
equations of Eq. (3. 6) are the three equations of motion. It is important to realize that 
we are not saying that all the forces actually apply at the center of mass. All that is 
claimed is that the motion of the center of mass can be determined by considering that 
all forces act there. 

3. 4 ANGULAR VELOCITY. 

3. 4. 1 We now discuss the fundamental concepts of angular motion. The angular velocity 

of the body axes about the translating- axes is the same as the angular velocity of the 

particles in the body rotating about some instantaneous axis. Since the angular velocity 

of a particle is a vector (paragraph 2. 10. 1) , the angular velocity of the body will also 

be a vector. The actual proof that it is a vector is beyond the scope of this report; 

however, we can intuitively see the truth of this fact. For example, if the body is 

moving such that its X- Z plane is rotating about its Y axis, the angular velocity 

vector lies along the Y axis. If the body is instantaneously rotating about any line, the 

instantaneous angular velocity will be along this line and will have direction according 

to the sense of the rotation. If applied forces were to cause the body to rotate about 

several axes of rotation simultaneously, the body angular velocity vectors about each 

line of rotation can be added vectorially into one total (resultant) angular velocity 

vector a; . Since any vector can be resolved into components with respect to any 

coordinate system, we can resolve the body angular velocity along the body axes. The 

components along the body axes X, Y , and Z are generally referred to as p, q, and r 

respectively; occasionally, to , co , w are used. Thus, 

x y z 

co - pi + qj + rk 
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the magnitude and direction of each force and the distance from the point of application 
to the center of mass. The rotation -causing force is related to a concept known as 
torque . The resulting rotation and associated angular velocity are involved in a con- 
cept known as angular momentum . A relation connecting the cause (torque) and the 
effect (angular momentum) called the moment equation is developed in subsequent 
paragraphs. First, however, it is necessary to discuss the concepts of torque and 
angular momentum. 
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The torque L about point P exerted by force F at point Q is given by the cross- 
product of distance PQ and force F . In other words , 




L = r x F 
P 



(3.7) 



where r = PQ . 



The magnitude of the torque is given bj' 
rF sin 6 =• Fd, where d is the perpendi- 
cular distance from P to the line in which 
Fi S- 34 ^he f orce ac ts. Torque can then be con- 

sidered as a force-moment. If P is fixed and both P and Q are on a planar body 
B of uniform thickness, then the force shown in the figure would cause the body to 
rotate about P in a counterclockwise direction. Thus, the torque vector would be 
parallel to and in the same direction as the angular velocity vector resulting from the 
application of only this one force. 
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The angular momentum H of a particle of mass m about point P is defined as the 
cross product of distance r from P to the particle and the linear momentum m f 
of the particle : 
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Since the angular momentum is the moment of the linear momentum about P , it is 
also frequently called the mo ment of momentum . (We define the moment of a vector 
| about a point as the cross product of the distance from the point to the origin of the 

J vector with the vector itself. In this case, one may consider the vector to be localized 

instantaneously along a specific line of action). 

The angular momentum of a rigid body about point P is the summation of the angular 
I momenta of all the particles of the body: 

J H p - f (?i x mjr.) (3.9) 

J 3.6 COMPONENTS OF ANGULAR MOMENTUM AND MOMENTS OF INERTIA. 

I 3.6.1 Before determining the instantaneous rotation of the rigid body caused by forces , 

we need to find expressions for the components along the body axes of the angular 

J : momentum of the body with respect to its center of mass. 

j Let I , j , and k be the orthonormal vectors of the body coordinate system (X, Y, Z) 

I which is rotating with angular velocity w = pi+qj+rk with respect to the 

t translating reference system (X Q , Y , Z ); the origin of both systems is the body 

! center of mass CM. Suppose P. is a particle of mass m. in the body at a distance 

r i = X i^ + y i^ + z i^ from the CM " The velocity of P. with respect to the transla- 
ting reference frame is given by Eq. (2. 12) as 



r. = (r.) u , + co x r. 
l x l'body i 

However, since the distance between P. and CM is always constant, P. has no 

i l 

velocity relative to the body axes , and the equation reduces to 
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Hence, H (the angular momentum of the body about CM) is given by 



H = 2 [ r. x m. (cj x r.) ] 
i l i l 



ft 



2 — - — - 

2 m. [r. co - (r. • w) r. 1 
I i i v l i J 



(3. 10) 






Let us now define a as the average density of the body . Then, the mass of the i^ 1 
particle can be written as 



m. = ctAV. 
i l 



(3.11) 



where A V is the volume element associated with the particle. Resolving H into 
components along the body axes, and, substituting Eq. (3. 11) into Eq. (3. 10), we have 



n 



n 



n 

r 



H 



x L 



r. p - (x.p + y.q + z.r)x. 



r. q - (x.p + y.q + z.r )y. 



2 jr. r - (x.p + y.q + Y^z. 



a AV. 
l 



o- A V. 



a AV. 
l 



P 
i ■ 

4.: 



2 2 2 2 

Grouping in terms of p, q, and r, and substituting x. + y. + z. for r. , we get 
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H 



p 2 (y 2 + z 2 ) a A V. - q 2 x. y. a A V. 
. w i i ' i H . l J i l 

l i 



-p 2 y. x. a A V. 
^ . l l l 

l 



-p 2 z. x. cr A V. 
.11 l 

l 



r 2 x. z. a AV. 
l l l 

l 



+ q 2 (x 2 + z 2 ) a A V. - r 2 y. z. a A V. 
n .i i' l .11 l 

l i 



q 2 z. y. cr A V. 
H . l 'l l 



+ r 2 (x 2 + y 2 ) ctAV. 
. v l J i ' l 



But the limit of 



2 (y 2 + z 2 ) a AV. 
l 



as the maximum of the volume elements A V. approaches zero is 

/ (y 2 + z 2 )crdV 
B 

(this symbol represents the triple integral over the three-dimensional region defined by 
the rigid body). Since all the sums above can be replaced by integrals, we have 



H 



2 2 

P / (y + z ) a dV - q / x y a dV 



B 



-p / y x a dV 
B 



B 



- r / x z a dV 
B 



2 2 

+ q / (x +z)crdV -r/yzadV 
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Replacing the integral expression 



I ■ 
t, 

n 



/ (y 2 + z 2 ) ffdV 
B 



by the symbol I , , replacing the integral expression 



/ x y a d V 
B 



by I , and replacing the remaining integrals in similar fashion, we get 
xy 



H = 



(pi -ql -rl ) 
v r xx H xy xz ' 

(-pi + q I -rl ) 

v ^ yx yy yz ' 



(-Pi 



zx 



q I + r I ) 
n zy zz 



H i + H j + H k 
x y J z 



(3.12) 
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Rewriting in matrix format, we have 



H = 



xx 



yx 



zx 



xy 



yy 



zy 



-i 



[ 

xz 




"p 




"H 

X 


[ 
yz 




q 


- 


H 

y 


zz 




r 




H 

z 



(3.13) 
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rt 



*xx ' *vv ' anc * * are ca -^ ec * ^ ie moments of inertia of the body about the body axes 
X , Y , and Z , respectively. The other I terms are called the products of inertia . 

The equations for the inertia terms are, of course, dependent on the geometric con- 
figuration of the body. The integral over the body with respect to the volume is replaced 
by the equivalent triple integrals written with respect to the body coordinate system.. 
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the origin (i. e. , the three mutually perpendicular axes about which all products of in- 
ertia vanish). That such directions exist is given by a theorem in classical mechanics 

I which will not be proved here. The moments of inertia taken about the principal axes 

are generally called the principa l mom ents of inertia at the point selected as origin of 

1 the axes. 

| 3. 7 MOMENT EQUATIONS FOP. BODY ORIGIN AT CENTER OF MASS. 

J 

3.7. 1 So far, we have obtained detailed expressions for angular momentum and have 
| shown that angular velocity is a component of angular momentum. The problem yet g 

to be resolved is that of finding the relation between the torques of impressed forces 
I and the resulting angular momentum. 

| Consider the equation of motion for each particle m. of a rigid body : 

i l 

9 m.r. - F. (3,34) 

1 1X1 

th - 

F. is the sum of the forces acting on the i particle with position vector r. from 
J x l 

l the origin O of some inertia! coordinate system. Taking the cross product of both 

sides of the equation with r. , we obtain 

| i 

J _ 

r. x m.r. '= r. x F. (3. 15) 

11111 

But note that 



—r (r. x m.r.) = i*. x m.r. + r. x m.r. - r. x m.r. 
dt v i ii i ii i ii i ii 

Substituting in Eq. (3. 15) , there results the fundamental equation 
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Summing over all particles of the body, 



thus 



Z I dt / „ ~ Z (3 Vo 



o i 



or 



dH 
o 



= L, 



o 

Equation (3. 17) stales that the time rate of change of the angular momentum of a rigid 
body about an incrtial origin O q is equal to the sum of all the torques impressed on 
the body. 

The equations for angular momentum and torque are much more tractable if these 
quantities can be expressed about the center of mass. The question then arises: Does 
the rate of change of angular momentum about the mass center equal the impressed 
torque about the mass center? The answer: Yes; regardless of the fact that the center 
of mass 5s moving with respect to inertia! axes (as v/e shall demonstrate). 

Once again the axes X , Y, and Z (Fig, 35) form a Cartesian coordinate system aUached 

to the rigid body B with origin taken as the center of mass CM . P. is a particle of 

mass m at distance R from the reertial origin O and at a distance r. from CM. 
1 i . o i 
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Fie. 35 



Let p be the position vector of CM with respect to O . Then, the angular momen- 
tum of B about O o [Eq. (3.9)] is 



H~ = 2 (R. x m.R.) 
j v i i i' 



2 Up + r ) x m (p + r,)] 

j 1 1 A 



2 (p x m.p) + 2 (p x m.r.) + 2 (r. x rn.p) + 2 (r. x m.r.) 
i l i lijii iiii 



p x (2 m.)p + p x (2 m.r.) + (2 m.r.) x p + 2 (r. x m.r.) 

i 1 ill ill i 1 1 J. 



(3.18) 
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But fl 

gives the position of the center of mass with respect to the body origin. Since the body H 

origin and the center of mass are coincident in this case, $-3 

2 m.r. = 
i i i 



Substituting these expressions into Eq. (3. 18), we obtain 



&*j 



and hence, ! 

t.; 

2 m.f . = |1 

i l x M 

Note, too, that I , 

2 (r. x m.r.) H 

1 ti 

is by definition the angular momentum II,,,,. of the body about its center of mass. » 



{' 



*..<•' 



H Q = p x mp + H CM (3.19) |j 

Differentiating Eq. (3.19), I : 

H Q = p x mp + p x mp + h T] 

o ' i \ 

• H Q - p x 2 F . + H (3.20) 
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since 



| tap - S F. 

< i x 

| bv Newton's second lav/ [Eq. (3.6)] , where 

j 

i 2 F. 

! i x 

is the sum of all the external, forces acting on B . The torque L„ exerted by these 

° 

! forces about O is 

• o 

L n = 2 (I. x F ) 

i u o i 1 * 

[ = 2 [(p + r.) x F. ] j 

. IIP j/ jj i 



2 (p x F.) + 2 (r. x F.) 



L ^ X f F i + L CM (3 " 21) 



where 



2 (r. x F.) 



is by definition the torque L about the center of mass. But, since we showed with 
Eq. (3. 17) that K Q = Lo • we can substitute Eq. (3. 20) into Eq. (3. 21); from which, 
we verify that 
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Thus, the rate of change of angular momentum of the body about the center of mass 
CM is equal to the sum of the torques about CM caused by the forces impressed on 
the body. For simplicity, we rewrite Eq. (3. 22) as 



n 
r: 



H - L 



If H is measured with respect to the body axes as in the preceding section, we have 



Li 



H = < H W + w x H 

where co is the angular velocity of the body with respect to the translating reference 
frame. Expanding by components, we obtain 



n 



H - H i + H j + H k + 
x y z 



i 



H 



J 

q 

H 



k 

r 
II 



(H + H q - H r) i + (II + H r - H p) j 
% x 7, y v y x z ' J 

+ (H + K p - H q)k 

v z y x' 



(3. 23) 



y 

n 



F: 
li 



Let 



L = Li + Lj+Lk 

x y J z 



where L , L , and L are the sums of the torques about bociv axes X , Y , and Z , 
x y z 

respectively. 
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since 



rap -- 2 F. 
• i i 



by Newton's second law [Eq. (3. 6) ] , where 



2 F. 

i i 

is the sum of all the external forces acting on B . The torque L r exerted by these 

u o 
forces about O is 
o 

L„ = 2 (1. x F.) 
°o i i i 



= 2 [(p * r ) x F ] 
l 11 

= 2 (o x F.) +2 (r. x F. ] 
i i> i i i' 



or 



L r . =■ p x 2 F. + L,„, {3, 21) 

°o i 1 CM 



where 



2 (f. x F.) 

is by definition the torque L„,, about the center of mass. But, since we showed with 
Eq. (3.17) that H = L Q . we can substitute Eq. (3. 20) into Eq. (3. 21); from which, 
we verify that 
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Thus , the rate of change of angular momentum of the body about the center of mass 
CM is equal to the sum of the torques about CM caused by the forces impressed on 
the body. For simplicity, we rewrite Eq. (3. 22) as 



r 

ii 

n 

Li 



H = L 



If H is measured with respect to the body axes as in the preceding section, we have 

i= Wbody^ xFl 

where co is the angular velocity of the body with respect to the translating reference 
frame. Expanding by components , we obtain 



H - 



H i + H j + H k + 

x y z 



j . 


k 


q 


r 


H 


H 



24 






li 

li 



n 
li 



Let 



(H + II q - II r) i + (H + H r - H p) j 
v x z^ y ' x y x z ' J 



+ ( H + H p - H q ) k 

' z y x ' 



(3.23) 



I; 






L = Ll+Lj+Lk 

x y J z 

where L , L , and L are the sums of the torques about bodv axes X , Y , and, Z , 
x y z ^ 

respectively. 
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I By equating components i , j , and k of L and H and substituting the expressions 

in Eq. (3. 12) for H , H , H , . . , , we obtain three equations: 
a x y 

L = p I - q I - rl + p I •- q I -rl + q ( - p I -q! + r I ) 

x xx "- xy xz r xx u xy xz ^ x *. zx H zy zz' 

I -r(-pl+ql- x I ) 

v c yx ■ yy - ? yz ' 

1 ...••• 

I L = - p I + q I -rl -pi +al -rl +r(pl -ql -rl) 

I y - yx ^ yy yz ^ yz " yy yz s ^ xx ^ xy xz ' 

1 -p(-pl-ql+rl) 

I ^ v r zx l zy zz ' 

i 

,, L = -pi - ci I + f I -pi -ql +rl +p(-pl + q I -rl ) 
I z * zx l zy zz ^ zx x zy zz * v ^ yx ^ yy yz ' 

I 

-q(pl -ql -rl ) 
i \ t- :KX ^ xy xz ' 

With angulax accelerations on. the left, we rev/rite the three equations (above) as 



pi -ql -rl = L -(pi -ql -rl )-q(-pI - ql +rl J 

1 xx H xy xz x x * xx H xy xz ' H s zx M zy zz ' 

+ r(-p I + q I - r I ) 

v * yx H yy yz ; 

pi + q I - r I = L - (p I + ql - r I ) - r(p I -ql -r I ) 

y x yy y z y y^ yy y z v * xx * xy xz 

4p(-pl -ql + r I ) 
F v *zx x zy zz ' 

•pi - ql + r I = L -(-pi -qi + r I ) - p ( -p I +ql -rl ) 

* zx H zy zz z l F zs H zy zz ' F v v yx H yy yz 

+ q(pl -ql -rl ) 
H v ^ xx H xy xz ' 

(3.25) 



f! 



3-80-61-2 

4 Jan 1963 



integrated to obtain the angular velocities p , q , and r . The derivatives of the direo- M 

tion cosines relating the body coordinate system to the translating reference frame can- 
then be computed as shown in paragraph 3.8.1. P 

If the body possesses certain geometric properties, the moment equations become n 

simplified. For example, if the body is symmetric with respect to all its axes, the prod- *--' 

uct of inertia terms and their derivatives vanish. An interesting analog to Newton's 
second law of motion results for a symmetric body of constant mass rotating with con- ; t 

stant velocity about only one of its axes. Let us say that this body is rotating with con- 
stant angular velocity p about its X axis;q and r are zero. Because the mass is ff 
1 constant, I is also zero. The moment equations reduce to one equation: 

If XX 

I p = L 
>oc x 

Note the similarity, in this case, to motion of the center of mass in the X direction: ^ 

n 

ma = F t i 

3. 8 CHANGE OF BODY ORIENTATION. \ I 

3. 8. 1 So far, we have obtained this relationship between the impressed forces and the 
angular motion of a body. We shall now show how the body orientation varies with the 
angular velocity. As stated In paragraph 3. 2. 1, the instantaneous orientation of the 
body with respect to the translating reference frame (X , Y , Z ) may be given by 
the orientation of a Cartesian coordinate system (X , Y , Z) to which the body is rigidly 
attached. In this case, the origin of the body axes X , Y , and Z is taken at the body 
center of mass. 







t! 

U 

I' 
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as the Euler angles. The orientation problem, then, is to determine how the transfor- 
mation varies with time as a result of the angular motion of the body. This can be 
determined by integrating the time rates of change of either the Euler angles or of the 
direction cosines; for reasons stated at the end of this section (3.8), it is the direction 
cosine approach which is developed. 



Before proceeding further, however, we need to expand upon the relation between the 

•f 

| time derivative of a vector A as measured in one' coordinate system, (which will be 



called the reference system) and the derivative of the same vector as measured in a 
coordinate system moving relative to the, reference frame. This relationship is devel- 
oped in paragraph 2. 11. 1 and stated in Eq. (2. 14), which we shall write in the more 
general form 

(#)..._-(£)_ ♦'♦ = ** < 3 - 26 > 



'reference relative 



where 



ldJ\\ 

•jT-j = the derivative of A as seen by an observer in the reference 



reference 



coordinate frame 



'm 



dt / 

relative 



= the derivative of A as seen by an observer in the frame 
moving relative to the reference frame 

p = the velocity of the origin of the moving (or "relative") frame 

co = the angular velocity of the moving frame relative to the refer- 

ence system. 

Note that, if the relative frame is only rotating and not translating, the velocity of its 
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t 



The reference system in Eqs. (3. 26) and (3. 27) is considered fixed in space only as 

seen by an observer in. the relative system. In point of fact, the reference system may 

be moving with respect to a third coordinate system, but Eq. (3. 27) is still a valid 

expression which relates one coordinate system that is rotating with respect to anotber. 

Suppose that the reference system is only translating with respect to an inertial system 

and that its axes remain parallel to the corresponding inertial axes: Then, the 

orientation of the relative system (X, Y, Z) with respect to the inertial system (X 1 , 

Y' , Z' ) is identical to the orientation of (X,Y, Z) with respect to the reference system 

(X„ , Y , Z) . ' In our previous discussion, we called frame (X , Y , Z ) the 
v ooo' r v o o o 

translating reference frame and actually stipulated that its origin coincide with 
the origin of the body axes X, Y, and Zi Under these conditions, Eq. (3.27) gives 
the derivative of A with respect to the translating reference system; adding the velocity 
of the common origin, we obtain the total derivative of A with respect to the inertial 
frame [Eq. (3, 26)] . 

We return now to the problem of finding the instantaneous orientation of the body axes 
with respect to the translating reference frame. Let i ,1 , k be the orthonormal 
vectors of the translating reference system, and let i , j , and k be the orthonormal 
vectors of the body-fixed coordinate system. As shown in paragraph 1.4. 1, these two 
sets of vectors are related by a matrix of direction cosines, which transforms coordin- 
ates in body axes to coordinates in reference axes : 



r-i 

i; 

v 

« . 

ii 



n 

i * 

Lj 



U 



n 



c c 

xx yx 

C C 

xy yy 



C C 

xz yz 



zx 




i 


zy 




3 


zz 




k 



(3.28) 
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Suppose that the body or "relative" coordinate system is rotating with angular velocity 
Zi =pi+qjH-rk with respect to the translating reference system (X o , Y o , 1^) . 
Then Eq. (3. 27) can be used to obtain the derivatives of \ , j Q , and k Q . Thus, 



r di 
_° 
, dt~ 



reference 



^7 



+ a) x i 



body 



W _L°\ 



dt / , , o 

/ body 



Therefore, 



= since I is fixed in magnitude 

and direction in the reference 
system. 



( 



dl \ 

_-°-J = i x w = 

/ body 



l 
C, 

P 



xx 



J 
C 

q 



yx 



k 
C. 



zx 



cr 



dt J 



body 



frC - a€ )1 + (nC - rC )j + (qC - pO )k (3,29) 
{ yx q zx ; u zx xx ; - 1 XH xx l yx ; x 



Sirnil arly , 
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anc 



./dk 



I 



dt 



body 



= (rC V9! - qC„)i + ( P C zz - rC xz )j + ( q C^ - pC_)k (3.31) 



yz x zz' 



xz *• yz 



Arranging Eqs. (3.29), (3.30), and (3.31) in matrix format, we have 



di 
c 

dt 

dj 

df 

d£ o 

dt 



o 



[rC yx- qC W fP C Z x~ rC xxJ ^ C xx-P C yx ; 



[rC yy" q V [pC W- r V [q C xy " p V 



[r C - q C ] [p C - r C ] [q C - p C ] 
I body • y 



But by straightforward differentiation of Eq. (3. 23), we obtain 



(3.32) 






Li 



di 

o 

dt" 

dj 
_£ 
dt 



dk 



dt 



body 



c c c 

xx yx zx 



c c c 

xy yy zy 



c c c 

xz yz zz 



k 



(3.33) 



Note that in differentiating Eq. (3. 28), the terms involving- 



m 



m , m 



r t i 
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Thus, by equating the matrices of Eqs. (3. 32) and {3. 33), we obtain a set of nine 

differential equations, which, when integrated, will furnish the direction cosines that 

give the orientation of the body coordinate system with respect to the reference frame. 

That is, for C and C (as an example), we obtain 
xx zy 

dC 



c 

XX 


dt 


= r C 

yx 


- q C 

n x> 


C 

zy 


dC 

Z Y 
dt 


= q C 

xy 


p C 

yy 



The direction cosine differential equations are thus functions of the direction cosines 
and the body angular velocity. They are continuous functions, whereas the Euler 



J angle differential equations (not derived here) are discontinuous and less tractable. It 

I 

also has been found a superior method (from a numerical computation viewpoint) to 

i obtain the transformation by integrating the direction cosine derivatives rather than by 

i 

I integrating tee Euler angle derivatives and then computing the transformation by the 

methods shown in Section 1.4. Errors obtained in the numerical integration of the 

| Euler angle rates may be compounded by roundoff errors and slight inaccuracies in the 

trigonometric calculations that are required to obtain the transformation from the 

j Euler angles. 

1 The direction cosine derivatives are smoother time functions than the Euler angle 

derivatives and, hence, can be numerically integrated with greater accuracy. The 

i properties of the direction cosine matrix \ the dot product of one row (or column) with 

a different row (or column) being zero, and the dot product of one row (or column) 
with itself being equal to one, as discussed in Section 1.4. I] can be used to check the 
accuracy of the integration or to reduce the number of equations being integrated. 
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3.9 MOMENT EQUATIONS FOR ARBITRARY BODY ORIGIN. 



3. 9, 1 It is important to remember that the moment equations derived in Section 3.7 
are valid only when torques are taken shout the center of mass. The following discus- 
sion develops the moment equations for the case in which torques are taken about a 
body origin which is a point of the body not at the center of mass. Once again, we 
employ the translating reference frame (X Q , Y , Z ) whose origin O coincides with 
that of the body axes X, Y, and Z . Let the position vector of the center of mass with 
respect to O be' r CM . P.. is an arbitrary point of the body at distance R. from O 
and distance f . from CM . 







V 

u 



u 



The angular momentum of the body about it 



f; 



r i 
% ■ 

n 






n 
I.; 

n 

Li 



u 



= 2 (R. x m R ) 
i 



n 
i, 



= V 



Sf(r CM+ r.)xm i( ? cM , ?.}] 



n 

i : 

i. J 
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But, 



2 m. r. = and 2 m. r = 
1 .11 ,.i 

I i i 

J by definition of the center of mass [Eq. (3.5)] , and 



I 



Z (r, x m. r.) 
.i l i' 

l 



I is the angular momentum H c , r of the body about the center of mass. Substituting 

these relations into Eq. (3.34) and solving for ll pfr , we obtain 



H CM "' H o r CM X m r CM 



.«> 



J Therefore , 

1 

} The torque of the body about O is 



r CM x Z F. + Z (r, x F.) 

i l 



or 



T = v v V 4 T 

o CM CM 



\ 
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Therefore, 

L CM = E o " f CM X f ( 3 - 36 ) 



r CM = < r CM>body + w Xr CM + ^ X < r CM>body + S X <° X ? CM> 



n 

f ! 



Since H CM •■= L CM » the right-hand sides of Eqs. (3.35) and (3.36) are equal, p 

giving 1 1 

5 o " r_ CM X 'CM = L o - ? CM X f < 3 " 37 ) ^ 

*"\ 

This expression [Eq. (3.37)] is the moment equation for the case in which forces {J 

and torques are measured about a body origin not at the center of mass: it is valid, 

f » 
of course, even when the center of mass is moving with respect to the body coordinate J j 

system . 

I; 

In order to write the detailed component moment equations along the body axes, we 

need to find an expression for r (the acceleration of the mass center with respect r 1 

to the translating reference frame centered" at O) . For this purpose, we use Eq. (2.15), '"' 

which becomes - _.. 

• r 



u 

(3.38) 

10 



Since the origins cf the body axes and the reference frame coincide, the acceleration ^ 

p [ Eq. (2, 15) ] of the body origin with respect to the translating reference frame is 
zero. The angular velocity of the body with respect to the reference frame is 



P 



oJ = p i + q j + r k |j 
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We now express each, of the terms in Eq. (3. 38) by components along the body axes: 



* r CM*body 



c J C J c 



(3.39) 



w X *CM = (qz e " ry e ):i + li-X c " pz e^ + (py c " qx e>k 



(3.40) 



2w X < ? CM>body = 2(qZ e "^V 1 + 2(rX c " pZ c }I + 2 ( py c " qX c )£ (3 " 41) 



w x (w x r CM ) 



(px c + qy c + rz c )p - o> x c 



(px c + qy c + rz c )q - w y c 



(px + qy + rz )r - w z 



Substituting these expressions into Eq. (3.38), we obtain 



(3.42) 



'CM 



x c + (q?J c " ry c ) + 2(qz c " ry c ) + (px e + qy e + rz c )p " w x c 

2 
y +(fx - p z ) + 2 ( r x - pz ) + (px + q y + r z ) q - co x 
J c c c c c 7 xr e lJ e c H c 

2 
z c + (Py c - qx c ) + 2(py c - qx c ) + (px Q + qy c + r z q ) r - w x c 



The second term on the right-hand side of Eq. (3. 37) becomes 



(3.43) 



r CM x F 



(Y F - z F ) i + (z F - x F ) 1 + (x F - y F )k 
' v c z c y ' x c x c z' - v ey •'ex 



(3.44) 



The three expressions in Eq. (3.25) represent the moment equations when H = L 
Ihus, we need to add the appropriate components of r x mr„ lf - r mr x F to 



riuff »/im 






the right-hand sides of Eq. (3, 25) in order to obtain Eq. (3. 37). This we will do to tr 
first expression in Eq. (3. 25), which represents the I component of Eq, (3.37): 

Pkx " qI xy " * J xz ~ L x " ^xx - qi xy ~ r W ~ <J<-P I ZX " ^y + ^ zz ) 

+ r(-pl yx+ q Tyy - rl yz ) +y c fm^ + (py c - q> c ) 

+ 2 (Py„ - qx ) + (px + qy + rz )r - o> 2 z ] 

V» W C» C/ Cr C 

' Z c [m ^c + ^ X e ~ P Z r> + 2 < r *n - Pi) + (px„ + qv 



r- 



«c-- 



PI 



+ rzjq - w 2 y 1 - (y F - z F ) « 

c' n •'c J w c z c y ; 

Collecting all p , q , and f terms and rearranging certain others, we obtain L 

P fl xx - m(y 2 c + z 2 c )] - qd^ - mx^) - f (1^ - mx^) (3.45) J 

c L x " < pi xx - ql xy " r W " q( ~ pI zx " qI xy + r W + T l~*\x + qI yy ' r V ■ 

t 

+ y c mfM c + 2(py c ' q V + r(px c + qy c + rz c )] ' z c m[y c + 2(rx c ' pi c } 

+ q(px + qy + rzj] -- (y F - z F ) 
c n * c c e z c y 

Performing the similar operations on the 2 nd and 3 rd expressions of Eq. (3. 25) we I' 

obtain the j and k components of Eq. (3.37): 

~ P(I xy " m Vc ) + q[1 yy " m(X c + Z? c )] ~ * (I yz " my c Z c } < 3 * 46 > 

- I, ■ - (-pi + q i _ r f ) - r (pi _ q j _ r j )+ p (- p i - ql + rl ) 

y yz yy yz' ^ xx H xy xz ; M ' zx * zy zz 

i- z c mfx c + 2(qz c - ry c ) + p(px e + qy c + rz c >] - x c m[z c + 2(py c - qx £ ) [_ 
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\ -P< T xz " m Vc> " *l\z - my cV +r[I z Z - m ^ x ; + y c )] < 3 - 47 ) 

1 = L z - ( - pi zx - qt zy + r W - P( - pI yx + qI yy " *W + q(pI xx " qI xy ~ r W 

i .... .... 

+ x c m[y c + 2 ( rx c - P z c ) + q (P x c + qy c + ri4 c )] " y c m[x c + 2(qz c " ry c ) 

i 

I + p{px + qy + rz )] - (x F - y F ) 

rv c 1 c c cy ex 

5 
I 

! Equations (3.45), (3.46), and (3.47) are the complete moment equations for the case 

in which forces and torques are measured about a body-axes origin not at the center of 

J mass. It is apparent that the moment and product-of- inertia terms on these equations 

i 

must be computed with respect to the body origin rather than the center of mass. No 

| general expressions can be written for the position, velocity, and acceleration of the 

1 mass center relative to the body origin (that is, for: x ,y , z ; x > , y , z ; and 

I x , y , i ) as these quantities are functions of the specific body configuration and 

I c c c 

I motion to be studied. 

J 3. 10 ACCELERATION OF A POINT FKED IN THE BODY. ' 

I 

i 3. 10. 1 In applications, it is frequently necessary to compute the acceleration of a body 

i _ 

point with respect to the center of mass, where the distance r (between point P and 

! the mass center CM.) remains constant. Equation (2.15), once again, gives the correct 

i _ _ 

expression except that, in this case, the velocity v and acceleration a of the point 
i with respect to the center of mass are zero, since the distance between P and CM is 

1 constant. Since the origins of the body axes and the reference axes coincide, the accel- 

eration p of the body origin with respect to the translating reference frame is zero. 
So, Eq. (2. 15) becomes 

r = co x r + w x (a; x r) 
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Lettinff 



w =pi+qj+rk and r = xi + yj + zk 



then 



r = 



qz - vy + (px + qy + rz ) p - (p 2 + q 2 + r 2 )x 



r x - pz + (px + qy + rz ) q - (p ' + q + r )y 

2 9 9 

Py - qy + (px + qy + rz)r - (p + q" + r )z 



or 



Fl 



n 

1 J 

n 

i.i 



r - 



9 2 

(~q" - r )x + (-f + pq)y + (q + pr)z 

2 2 

(r + pq)x + (-r - p )y + (-p + qr)z 

2 2 

(-q + pr)x + (p + qr)y + (-p - q )z 



(3. 48) 



f •» 



To obtain the total acceleration of P , one must add to this the acceleration of the 
center of mass with respect to the inertial coordinate system being used. 



3. 11 EULER ANGLE DERIVATIVES. 



n 



3. 11. 1 As discussed in Section 3.8, the time variant transformation matrix of rotation 
between two coordinate systems can be obtained by numerically integrating the time 
rates of change of the direction cosines. It is also possible to integrate the time deriv- 
atives of the Euler angles and then obtain the transformation matrix directly from the 



f] 
U 

f| 
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Let (X , Y , Z Q ) be a reference frame and let (X , Y , Z) be another reference 
frame obtained from (X o , Y Q , Z ) by three successive rotations. For this example, 
let us take the displacements to be a roll <p , followed by a pitch , and, finally, by 
yaw if> . Associated with each angular displacement is an angular velocity <p , 6 , 
and 4> , respectively. If a rigid body is affixed to the coordinate system (X , Y , Z), 
the body will have an angular velocity 53 equal to the vectorial sum of the three sim- 
ultaneous velocities , 8 , and ip ; i.e. , 



(0 = 0+04-$ 



fsP 



V J 



However, the body angular velocity is normally obtained by integrating the moment 
equations to obtain p , q , and r (the components of 5 along each of the X , Y , Z 
axes). Thus, we have the relation 



o) =pi+qj+rk=0 + 0+$ 

By taking the components of <f> , 6 , and $ along the body axes, we have expressions 
for obtaining <p , 6 , and ip in terms of p , q , r , <p , , and 4' • -Let us now ex- 
amine the three successive rotations as illustrated in Fig. 37. 




Y 




_ *». y« , v 




-^* y 



/Ml 
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The angular velocity <p (due to rotation <t>) is normal to the Y -Z plane rnd falls 
along X q . The components of in coordinate system (X , Y , Z) are obtained by 
applying the transformation from frame (X , Y , Z ) to frame (X , Y , Z). Thus, 



* ; 
U 



*X, Y, Z 

coords . 



cos 4> sin ip ' 



Jin $ cos ^ 



cos 6 -sin 



sin 6 cos £ 



cos <p sin <jt 



<P 



-sin (•) cos 



d> 



f! 



cos if' cos 6' 
sin ;// cos 6 
cp sin 



(3.49) 



n 



The angular velocity (due to rotation 6) is normal to the X* - Z' plane and falls 
along Y' . Tho components of in system (X , Y , Z) are obtained by trans for mati 
from coordinates X' , Y' , and Z' : 



ion 



L J 



X , Y , Z coords. 



cos i[> sin i[> 



-sin ip cos i/ 1 



cos 6 



sin f 



sin 6 




"o" 










cos 




_0_ 



: 



["$ sin til 



(3.50) 



fi 
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The angular velocity i (due to rotation ty ) is normal to the X"-Y M plane and falls 
along Z" . The components of J in system (X , Y , Z) are obtained from coordinates 
X" , Y" , and Z" . 



#■ 



X , Y , Z coords. 



cos ip sin 4' 



-sin 4' cos # 





"o" 












rr 







i. 




A 



(3.51) 



Finally, the total angular velocity u> is obtained as the sum of Eqs. (3.49), (3.50), 
and (3.51): 



w = 



4> cos ^ cos + 6 sin # 



- $ sin # cos 0+0 cos # 



<f> sin 



+ 4> 



(3.52) 



This system, of equations can be solved simultaneously for $ , 8 , and 4> » yielding 



<P 



p cos J) - q sin j& 

cos 



- p sin ^ + q cos $ 



(3.53) 



ib - r - <b Bin 



Numerical integration of these time rates of change will give the Euler angles for rota- 



TT',- ,™ ,1 J .eC,- 



^r -/vi^fi^* ,4i.piV 



